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Determining the quantum circuit complexity of a unitary operation is an important problem in quantum 
computation. By using the mathematical techniques of Riemannian geometry, we investigate the efficient 
quantum circuits in quantum computation with n qutrits. We show that the optimal quantum circuits are 
essentially equivalent to the shortest path between two points in a certain curved geometry of SU(3"). As an 
example, three-qutrit systems are investigated in detail. 

Due to the quantum parallelism, quantum computers can solve efficiently problems that are considered 
intractable on classical computers 1 , e.g., algorithm for finding the prime factors of an integer 2-3 and 
quantum searching algorithm 4 . A quantum computation can be described as a sequence of quantum 
gates, which determines a unitary evolution U performed by the computer. An algorithm is said to be efficient if 
the number of gates required grows only polynomially with the size of the problem. A central problem of quantum 
computation is to find efficient quantum circuits to synthesize desired unitary operation 17 used in such quantum 
algorithms. 

A geometric approach to investigate such quantum circuit complexity for qubit systems has been developed 
in 5 7 . It is shown that the quantum circuit complexity of a unitary operation is closely related to the problem of 
finding minimal length paths in a particular curved geometry. The main idea is to introduce a Riemannian metric 
on the space of fi-qubit unitary operations, chosen in such a way that the metric distance d(I, U) between the 
identity operation and a desired unitary 17 is equivalent to the number of quantum gates required to synthesize 17 
under certain constraints. Hence the distance d(I, 17) is a good measure of the difficulty of synthesizing 17. 

In fact, d- dimensional quantum states (qudits) could be more efficient than qubits in quantum information 
processing such as key distribution in the presence of several eavesdroppers. They offer advantages such as increased 
security in a range of quantum information protocols 812 , greater channel capacity for quantum communication 13 , 
novel fundamental tests of quantum mechanics 14 , and more efficient quantum gates 15 . In particular, hybrid qubit- 
qutrit system has been extensively studied and already experimentally realized 1617 . The higher dimensional version 
of qubits provides deeper insights in the nature of quantum correlations and can be accessed by encoding qudits in 
the frequency modes of photon pairs produced by continuous parametric down-conversion. 

In particular, the three-dimensional quantum states, qutrits are of special significance. For instance, in the 
state-independent experimental tests of quantum contextuality, three ground states of the trapped 171 Yb + ion are 
mapped to a qutrit system and quantum operations are carried out by applying microwaves resonant to the qutrit 
transition frequencies 18 . The solid-state system, nitrogen-vacancy center in diamond, can be also served as a qutrit 
system, in which the electronic spin can be individually addressed, optically polarized, manipulated and measured 
with optical and microwave excitation. Due to its long coherence time, it is one of the most promising solid state 
systems as quantum information processors. 

In this paper we study the quantum information processing on qutrit systems. We generalize the results for 
qubit-systems 7 to qutrit ones. The efficient quantum circuits in quantum computation with n qutrits are inves- 
tigated in terms of the geometry of SU(3"). Three-qutrit systems are investigated in detail. Compared with the 
results for qubit systems 7 , our results are more fined, in the sense that by using enough one- and two-qutrit gates it 
is possible to synthesize a unitary operation with sufficient accuracy. While from 7 , it is not guaranteed that the 
error of the approximation would be arbitrary small. 



Results 

A quantum gate on n-qutrit states is a unitary operator 17 i 
H(t) according to the Schrodinger equation, 



SL7(3") determined by time-dependent Hamiltonian 
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[7(0) =1, U(T) = U. 



(1) 



For qutrit case the Hamiltonian H can be expanded in terms of the 
Gell-Mann matrices. As the algebra related to the n-qutrit space has 
rather different properties from the qubits case in which the evolved 
Pauli matrices have very nice algebraic relations, we first present 
some needed results about the algebra su(3"). 

Let A,-, i = 1, 8, denote the Gell-Mann matrices, 
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Let 



a?=i®--® A*®---®/ 



be an operator acting on the ath qutrit with A k and the rest qutrits 
with identity I. The basis of su(3") is constituted by {A s }, s = 1, n, 
where 



H(« 



1 < < a 2 < ... < a, s n, 1 < ^ < 8. A s stands for all operators 
acting on s qutrits at sites ai, a 2 , . . ., a s with Gell-Mann matrices A^ , 
Afc 2 , . . . , Ai s respectively, and the rest with identity. We call an element 
in {A s } an s-body one. By using the commutation relations among 
the Gell-Mann matrices, it is not difficult to prove the following 
conclusion: 

Lemma 1 All s-body items (5 a 3) in the basis of su(3") can be 
generated by the Lie bracket products of 1-body and 2-body items. 
In the following the operator norm of an operator A will be defined 

by 



||A|| = max \\Ax\\ 

11*11 — 1 



(2) 



which is equivalent to the operator norm given by <A, B> = trA^B. 
The norm of above Gell-Mann matrices satisfies ||A,|| = 1, i = 1, 
2 

7, and || A 8 || = —=. If we replace A 8 with — Ag, the Gell-Mann mat- 

V 3 2 
rices are then normal with respect to the definition of the operator 
norm, and the basis of su(3"), still denoted by {A s }, is normalized. 

A general unitary operator U e SU(3") on n-qutrit states can be 
expressed as U = l/ 2 . • • 14 f° r some integer k. According to Lemma 
1, every 17, acts non-trivially only on one or two vector components 
of a quantum state vector, corresponding to a Hamiltonian Hi con- 
taining only one and two-body items in {AJ, s = 1, 2. 



The time-dependent Hamiltonian H(t) can be expressed as 

G G 

where: (1) in the first sum X,/, a ranges over all possible one and two- 
body interactions; (2) in the second sum X", a ranges over all other 
more-body interactions; (3) the h a are real coefficients. We define the 
measure of the cost of applying a particular Hamiltonian in synthes- 
izing a desired unitary operation 17, similar to the qubit case, 

(3) 



F(H) = JJ2 h l+P 2 Y, h2 °> 



where p is the penalty paid for applying three- and more-body items. 

Eq. (3) gives rise to a natural notion of distance in the space SL/(3") 
of n -qutrit unitary operators with unit determinant. A curve [17] 
between the identity operation I and the desired operation U is a 
smooth function, 

't/:[0,fr]-SU(3») 
[/(0)=/and U(t f ) = U. 

The length of this curve is given by d{{ U] ) = $ dtF{H{ t) ) . As d( [U] ) 
is invariant with respect to different parameterizations of [17], one 
can always set F(H(t)) = 1 by rescaling H(t), and hence 17 is generated 



defined by 



at the time tf = d([U]). The distance d(I, U) between I and L7 is 
d(I,U)= mind\U}. (5) 

V[l/] 



The function F(H) can be thought of as the norm associated to a 
right invariant Riemannian metric whose metric tensor g has 
components: 

{0, if (T^T 
1, if (t^t and a, x is one- or two-body . (6) 
p 2 , if a x and a, x is three- or more-body 

These components are written with respect to a basis for local tangent 
space corresponding to the coefficients h a . The distance d(I, U) is 
equal to the minimal length solution to the geodesic equation, {dHI 
dt, J) = i(H, [H, J]). Here (\ •) is the inner product on the tangent 
space su(3") defined by the above metric components, and / is an 
arbitrary operator in su(3"). 

From Lemma 1 in the basis { AJ of sm(3"), all the q-body items (q S 
3) can be generated by Lie bracket products of 1-body and 2-body 
items. To find the minimal length solution to the geodesic equation, 
it is reasonable to choose such metric (6), because the influence of 
there- and more-body items will be ignorable for sufficiently large p. 
It is the one- and two-body items that mainly contribute to the 
minimal geodesic. 

We first project the Hamiltonian H(t) onto H P (t) which contains 
only one- and two-qutrit items. By choosing the penalty p large 
enough we can ensure that the error in this approximation is small. 
We then divide the evolution according to H P (t) into small time 
intervals and approximate with a constant mean Hamiltonian over 
each interval. We approximate evolution according to the constant 
mean Hamiltonian over each interval by a sequence of one- and two- 
qutrit quantum gates. We show that the total errors introduced by 
these approximations can be made arbitrarily smaller than any 
desired constant. 

Let M be a connected manifold and T> a connection on a principal 
G-bundle. The Chow's theorem 19 says that the tangent space M q at 
any point q e M can be divided into two parts, the horizontal space 
H q M and the vertical space V q M, where M q = H q M © V q M and 

V q M = Q ( 0 denotes the Lie algebra of G). Let {A*} be a local frame 
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x?,xt 



of H q M. Then any two points on M can be joined by a horizontal 
curve if the iterated Lie brackets X?, \X^ , ■ 

[ tk [ 'it- 1 ' 

evaluated at q span the tangent space M q . 



Lemma 2 Let p be the penalty paid for applying three- and more- 
body items. If one chooses p to be sufficiently large, the distance d{I, 
U) always has a supremum which is independent of p. 

Proof. As SU(3") is a connected and complete manifold, the tangent 
space at the identity element I can be looked upon as the Lie algebra 
su{3"). For a given right invariant Riemannian metric (6), there exists 
a unique geodesic joining / and some point U e SU(3"). With the 
increase ofp, the distance d(I, 17, p) of the geodesic joining J and U e 
SU(3") increases monotonically. 

On the other hand, according to Lemma 1, 1-body and 2-body 
items in the basis {A s } can span the whole space su(3") in terms of the 
Lie bracket iterations. Under the metric Eq.(6), from the Chow's 
theorem we have that the horizontal curve joining J and 17 e 
S17(3") is unique, since the subspace spanned by 1-body and 2-body 
items is invariable. Or there exists such a geodesic that its initial 
tangent vector lies in the subspace spanned by 1-body and 2-body 
items. Hence the distance d(I, U,p) has a sup d 0 which is independent 
ofp. 

Lemma 3 Let H P (i) be the projected Hamiltonian containing only 
one- and two-body items, obtained from a Hamiltonian H(t) gen- 
erating a unitary operator U, and Up the corresponding unitary 
operator generated by H P (t). Then 



\U-U P \\< 



3"d([U}) 



(7) 



where j | • j j is the operator norm defined by (2), and p is the penalty 
parameter in (6). 

Proof. Let U and V be unitary operators generated by the time- 
dependent Hamiltonians H(t) and /(f) respectively, 

dU dV 
— = -MU, — = -iJV. 
dt dt 

By integrating above two equations in the interval [0, T], we have 



U-V = 



i{JV-HU)dt, 



where U{T) = 17, V(T) = V and 17(0) = V(0) = I have been taken 
into account. 
Since 



dVU 
dt 



= (-iJV)*U+V*(- iHU) =iV*[J-H)U, 



we have 



V*U-I= -i 



V*(H-J)Udt. 



Using the triangle inequality and the unitarity of the operator 
norm || • ||, we obtain: 

rT 



\U-V\\ = || V*(17- V)|| = || V* L/-I|| = 



V*{H-J)Udt 



dt\\V*(H-J)U\\ = 



dt\\(H- 



\\H\\ = 



<X>|<3yfc?+fcf+---+/^ = 3lJf|| 



Moreover, if H contains only three- and more-body items, we have 



F(H) = J P 2J2K=P\\H\\ 



Therefore 

d([U]) = 



dtF(H(t)) 



> 



dtF(H{t)-H P {t)) = 



pdt\\H(t)-H P (t) 



- 3 « 



dt\\H(t)-H P (t)\\>^\\U-U P \\, 



which gives rise to (7). 

Remark. From Lemma 3, by choosingp sufficiently large, sayp = 9", 
we can ensure that ||U — Up|| £ d([U])/3". Moreover, since the 
distance d(I, 17) is defined by d{I, 17) = min^LrjcftU], Lemma 3 also 
implies that \\U - U P \\ < d(I, U)/3". 

Lemma 4 If U is an n-qutrit unitary operator generated by H(f) 
satisfying ||H(f)|| £ c in a time interval [0, A], then 

|| U- exp(-iH)A||<2(e cA -l-cA)=0(c 2 A 2 ), 



where H e 



1 



dtH(t) is the mean Hamiltonian. 



Proof. Recall the Dyson series 20 : 

00 



f A 


■h 


dh 


dt 2 --- 


Jo 


0 



dt m H(h)H(t 1 )---H(t m ). 



We choose f, £ A/(z + 1) and set the first term in the above series 
to be I. Hence the second term is ( — i) J 0 H(ti)dti = — iAH. We 
have 



, ( — (HA) 2 (-iHA) m 

I+( — iHA) + — -+ ■■■ + y - ^ + 

2 m\ 



E 

m = 2 

CO 

^E 



( — iHA) 



( — (HA) 



CO ,A rf, rtn-i 

V(-i) m dtA dt 2 - \ dt m H(h)H 
m = 0 Jo Jo Jo 

co eA eti [t m — i 

y)(-i) m * «&••• dt, 

m= 2 Jo Jo Jo 



H(h)H(t 2 ) ■ ■ ■ H(t,„) 



"A m c m A" 
ml m! 



+ f dh V dh-- - [" 1 dt m \\H{t 1 )H{t 2 )---H{t m )\\ 
Jo Jo Jo 

2(e cA -l-cA), 



where we have used the standard norm inequality | |XT| | £ | |X| | 1 1 Y\ \ , 
the condition ||H(f)|| £ c, ^ dt x J 0 '' dt 2 ■ ■ ■ J 0 '"' _1 dt m = A m /ml and 

A A A 



A . Am 

-...-A=— . 

mm— I 2 m\ 



Proposition 1 If A and B are two unitary operators, then 



The Euclidean norm of the Hamiltonian H = h a a is given by 
|H|| 2 = h}- From the Cauchy-Schwarz inequality, we have 1 Proof. We begin with N = 2. It is easy to verify that 
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||A 2 -B 2 || = \\A 2 -AB + AB + B 2 \\ < ||A(A-B)|| + 
||(A-B)B||=2||A-B|| 

Now suppose that this inequality holds for N — 1, N a 3, i.e., 

| | A n-i _ b n-i|| < (N - - B\\. Then for N we have 



< ||A N - 1 (A-B)|| + || {A N - 1 -B N - 1 )B\\ 

= ||A-£|| + ||(A N - 1 -B N - 1 )|| 

<||A-B|| + (N-1)||A-B|| =N||A-B||. 

Lemma 5 Suppose H is an w-qutrit one- and two-body Hamiltonian 
whose coefficients satisfy \h a \ £ 1. Then there is a unitary operator 
U A which satisfies 

II „-iHA tt || ^ „ „2 a3 
e — Ua < Q « A , 

and can be synthesized by using at most c^/A one- and two-tribit 
gates, where q and c 2 are constants. 

Proof. We need a modified version of the Trotter formula 1 : let A and 
B be Hermitian operators, then ^ A+B ^ = e iAA 'e ,BA ' + 0(Af 2 ). We 
divide the interval [0, A] into N = 1/A intervals of size A 2 . In every 
interval, we define a unitary operator 

U A2 = g—ihmit? e —ihiait? . _ . e -ih L o L A 2 

There are L = 32m 2 — 24n = 0(n 2 ) one- and two-body items in H. 
From the modified Trotter formula, there exists a constant c 2 such 
that 



-iHA 2 



i(h l oi+h 2 <j 2 -t \-h L o L )A 2 



-ih L a L A l 



<c 2 n 2 A\ 
By using Proposition 1, we have 

j | e - iHA _ v n j | < N j | e -,HA 2 - l/ A2 ] | < C2 N« 2 A 4 = c 2 « 2 A 3 . 

It means that one can approximate e~' H& by using at most Ncin 2 = 
qw 2 /A quantum gates for some constant C\. 
From the above we have our main result: 

Theorem 1 Using 0(n K d(I, C7) 3 ) (KeZ) one- and two-qutrit gates it is 
possible to synthesize a unitary U A satisfying 1 1 U — U A \ \ £ c, where c 
is any constant. 

Theorem 1 shows that the optimal way of generating a unitary 
operator in SU(3") is to go along the minimal geodesic curve con- 
necting I and U. As an detailed example, we study the three-qutrit 
systems. In this case the right invariant Riemannian metric (6) turns 
out to be a more general one 21 , 



tr(HQ(J)) 
2x3 2 ' 



where Q(J) =sS(J) +T(J) +pQ (J),p is the penalty parameter and s 
is the parameter meaning that one-body Hamiltonians may be 
applied for free when it is very small, S maps the three-qutrit 
Hamiltonian to the subspace containing only one-body items, T to 
the subspace containing only two-body items, and Q to the subspace 
containing only three-body items. According to the properties of the 
Gell-Mann matrices, they satisfy [S,T] c T, [S,Q] c Q, [T,Q] <=T. 

Set L = Q(H), S = S{L), T = T(L) and Q=Q{L). From the geo- 
desic equation L = i[L,J-[L)], where J-=Q~ 1 , we have 



( S = 0 

t=i 



s+l i--)q,t 

[S,Q], 



(8) 



which gives rise to the solution 

fS(t) = S 0 
T(t) = e''^ 1 - s_1 ) s » e 't(i-p-')(s 0 + Qo) 



T 0 e 



')s„ 



(9) 



s Q(f) = e «(P" 1 -' i - 1 )SoQ oe -«(p- 1 -s- 1 )So. 



where S(0) = S 0 , T(0) = T 0 and Q(0) = Q 0 . 

The corresponding Hamiltonian H = Q~ l (L) has the form: 

H(t) = -S(t) + T(t) + - Q(t). According to the assumption (H(t), 

S ? tr(S 2 ) tr(T 2 ) 
H(t)) = 1 for all time f, we have r <s, r- < 1, and 

rr(Q 2 ) 1 2X3 2X3 

< p. The term - Off) in H(f) is of order p 1/2 and hence can 

2 x 3 2 p 

be neglected in the largep limit, with an error of order tp~ m . Also the 
term containing^ -1 in the exponentials of Tcan be neglected with an 
error at most of order i 2 (s 1,2 p~ 1 + p~ 1/2 ). Therefore one can define an 
approximate Hamiltonian 



H(t) = -S 0 - 



j-ifc-'SogiftSo + QoW e -'t(S 0 + Q 0 ) e its- l S 0 



The corresponding solution U(t) of the Schrddinger equation 
satisfies 



U(t) 



-U(t)\\<0{tp- l l 2 + t 2 {s 1 ' 2 p- 1 +p- 1 ' 2 )). 



Denote V = e -'' f (So+Qo) e i<s % fj. Then y=-!(S 0 + T 0 
and V = e " '' f(s ° + T ° + Qo) . Thus we have 



-Qo)V 



U(t)' 



:e -its- 1 So e it(S 0 + Q a ) e -it(S 0 + T a + Q 0 ) _ 



Generally one can expect that S 0 + Qo is much lager than T 0 , and 
So + Qo is nondegenerate. U can be simplified at the first-order 
perturbation, 



U(t)- 



(To) 



where 1Z, 



(So + Qo) 



To) denotes the diagonal matrix by removing all the 
off-diagonal entries from T 0 in the eigenbasis of S 0 + Qo- Therefore 
we see that it is possible to synthesize a unitary (7 satisfying 
|| U(t) — 17(f) || < c, where c is any constant, say c = 1/10. 

Discussion 

We have investigated the efficient quantum circuits in quantum 
computation with n qutrits in terms of Riemannian geometry. We 
have shown that the optimal quantum circuits are essentially equi- 
valent to the shortest path between two points in a certain curved 
geometry of SU(3"), similar to the qubit case where the geodesic in 
S!7(2") is involved 7 . As an example, three-qutrit systems have been 
investigated in detail. Some algebraic derivations involved for qutrit 
systems are rather different from the ones in qubit systems. In par- 
ticular, we used (2) as the norm of operators. The operator norm of M 
used in 7 is defined by | |JVf| | ! = max^j,^ =1 {|(i/'|M|i/')|}, which is not 
unitary invariant in the sense that ||M| |i = 1 1 UM\ | x = | \M (7]|i is not 
always true for any unitary operator U. For instance, consider 



M = 



0 1 

0 0 



and 17 = 



One has \\M\li = 1/2. 



However, ||M = 1/2 + V2/4. Generally, from Cauchy- 
Schwarz inequality one has | |M| | j < | |M| | . If AitM = I or hf< = M, 
then HiVflU = 
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Moreover, the final results we obtained are finer than the ones in 7 . 
Our result shows that if k in formula ^ =n d(I,U) is taken to be 

sufficiently large, 1 1 U — U A \ \ can be sufficiently small. However, the 
approximation error estimation in 7 reads 



|| [ /-L/,||<^ + 2^( e (V^_( 1 
+ c 2 d{I,U)n 4 A 2 . 



First, since d(I, U) is dependent of the penalty parameter p, there 
should exist a p-independent bound to guarantee that 2"d(I, U)/p 
is small for sufficiently large p. Second, if one chooses A as scale 1/ 
n 2 d(I, 17), the sum of the last two terms of the right hand side is 9/2 + 
c 2 /d(I, U) + O. Therefore the scale should be smaller, for example, 1/ 
n k d(J, U) and k > 3. As A takes the scale of l/n 2 d(I, U) in 7 , it can not 
guarantee that the error in the approximation could be arbitrary 
small. 

Due to the special properties of the Pauli matrices involved in 
qubit systems, many derivations for qubit systems are different from 
the ones for qutrit systems. Nevertheless, the derivations for qutrit 
systems in this paper can be generalized to general high dimensional 
qudit systems. 

Methods 

In deriving Theorem 1, we use Lemmas 2-5. Let H(t) be the time-dependent nor- 
malized Hamiltonian generating the minimal geodesic of length d(I, U). Let H P (t) be 
the projected Hamiltonian which contains only the one- and two-body items in H(t) 
and generates Up. According to Lemma 3, they satisfy 



\\U-U P \\< 



i"d(I,U) 



(10) 



Divide the time interval [0, d(I, U)] into Nparts with each of length A = d(I, U)/N. Let 
U J P be the unitary operator generated by H P (t) in the jth time interval, and U J M be the 

_ l «a 

unitary operator generated by the mean Hamiltonian H— — dtH P (t). Then 

using Lemma 4 and inequality H p (f) <4\/2n we have 



||l/i-[7i i |<2(e 4 ^" A -(l + 4\/2«A)). (11) 
As F(H) is scaled to be one, F(H) = ^E,- K +P 2 E„" K = L one has 



F(Hp(t)) = VT7K<^ Hence 



l|Hp(f)|| = 



<Y^\hA<L\]h\ + hl 



-flf <4V2«, 



where L = 32m 2 - 24« is the number of one- and two-body items in H(t), i.e. the 
number of the terms in H P (t)). 

Applying Lemma 5 to IV on every time interval, we have that there exists a unitary 
U' A which can be synthesized by using at most C]_n 2 /A one- and two-qutrit gates, and 
satisfies 



\uL-u'\ 



-'ffA_^||< C2 „2 A 3 



(12) 



Up and U A can be generated in terms of U p and U 3 A , respectively. We show how to 
generate Up by use of ^ below. First, U p can be generated by Hp: 



dU p 
dt 



-iH p (t)Up(t), U p (0)=I 

with Up(A) — Up. The unitary operator Up generated by Hp satisfies 
dU 2 p (t) 



dt 

which can be transformed into 



-iH 2 p(t)U p (t), C£(A) = [& 



dU 2 p Up(t) 
dt 



= -iHp(t)U 2 (t)U P (t), Up-U P {0) = U P , 



U P =U"U" 



can be generated similarly. 



• Up generated by the Hamiltonians H l p (t), H 2 p (t), ■ ■ ■ ,Hp (t).U A 



Therefore 

\\Up-U a \\ 

= ||«-' 



■ ui-uyu* 



■u 



<\\UpU 

= \\u 



rWrrN-1 

lp Up 



■ u\ -c/^- 1 • • • u A +u» u?- 1 - ■ ■ ul-u» u"- 1 ■ ■ ■ u\\ 

■ ^(i/i-i/Dll+IKi^D*- 1 • • • vi-vw-i -. U\)V A \ 



r l|l_ L ll rr N rr »-l... [/ 2_ [/ M [/ N-l. 



\ + \\ul-u\\\ + ---+\\u%-u* 



^\\ U P- U A„ 

From (10), (11) and (12) we obtain: 

||!/-[/a||<||1/-[/p|| + ||[/p-[/a|| 

3"d(I,U) Am , , I 



••Pill- 

N 



p 

. rd(i,u) 
p 

. 3"d(I,U) 



+ 2N(e 4V2 " & -(l 



+ U'-U' 



4V2«A^ +c 2 Nn 2 



(13) 



= +2 d ^ J l (e^- (l +4vW)) +c 2 d(I,U)n 2 A 2 

= 3 " d V' U ) +2c 0 d(I,U)n 2 A + c 2 d(I,U)n 2 A 2 , 

where e is/2nL - [\ + i-/2nA^ = 0(« 2 A 2 ) and c 0 is a constant. 

As mentioned in Lemma 2, the distance d(I, U) has a sup d 0 for sufficiently large p. 
For example, we choose a suitable penalty p so that d(I, U, p) satisfies 

8d 1 

<d(I,U,p) <d 0 . If we choose A to be sufficiently small, e.g. — —n k d(I,U) with k 

sufficiently large, \ \U - U A \\ will be sufficiently small, 



\\u-u A \\< 



< 



Td(I,U) 
P 

3"do 



-2c Q n 



-2c„«-( t - 2 '- 



c 2 n 



-(2k-2) 



-n 



d(I,U) 

-(2k-2) 



(14) 



As we need c^/A one- and two-body gates to synthesize every U } A , we ultimately 
need N— — — — ^ — - — cid(I,U) 3 n 2k+2 one- and two-body gates. 
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